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Abstract. We present explicit generators D\, . . . , _D„ of an algebra of commuting differ- 
ence operators in n variables with trigonometric coefficients. The algebra depends, apart 
from two scale factors, on five parameters. The operators are simultaneously diagonal- 
ized by Koornwinder's multivariable generalization of the Askey- Wilson polynomials. For 
special values of the parameters and via limit transitions, one obtains difference operators 
for the Macdonald polynomials that are associated with (admissible pairs of) the classical 
root systems: An~i, Bn, Cn, D„ and BCn- By sending the step size of the differences to 
zero, the difi'erence operators reduce to known hypergeometric differential operators. This 
■^.^ limit corresponds to sending q ^ 1; the eigenfunctions reduce to the multivariable Jacobi 



On 



polynomials of Heckman and Opdam. Physically the algebra can be interpreted as an 



integrable quantum system that generalizes the (trigonometric) Calogero-Moser systems 
related to classical root systems. 
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1 Introduction 

Over the past few years, progress has been i 
in more than one variable. It has turned 
orthogonal polynomials depending on only < 
ables. Such generalizations can be viewed 
root system R, there exist associated fam 
ber of variables coincides with the rank n 
Jacobi polynomials related to root system: 
I HO, |H1|. Recently, a more elementary a( 
in [|H2|1. For R = BCi the Heckman-Opd 
Jacobi polynomials in one variable. In a y( 
troduced q-versions of th e Heckman-Opda: 
of these results and [ [M3| ] for lectures devo 
then Macdonald's polynomials coincide wi 
information on continuous q- Jacobi polync 
results pertaining to the root system BCn h 
finds BCn-type multivariable versions of th 
one- variable case is recovered by specializin 
A crucial ingredient in the constructior 
operator of which the polynomials are eig( 
Jacobi polynomials this operator is a sect 
which is named hypergeometric differential 
polynomials the relevant operator is an an. 
transformation turns the hypergeometric I 
adjoint with respect to Lebesgue measure, 
upon the latter PDO as being the Hamih 
dimension one. Such quantum systems have 
literature; they are known as (generalized] 



Su; 
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The systems studied in and Sul 
the survey paper [3P2 for the generalizat 
the corresponding classical systems has al 
(arbitrary R). 

It has been shown by Heckman and Opd 
tor is but one member of an algebra of com 
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polynomials as their joint eigenfunctions |H0, HI, 01, 02 1 
independent PDO's 



This algebra is generated by n 



This state of affairs can be expressed by saying that the correspond- 
ing (generalized) Calogero-Moser system is quantum integrable. For R = An-i, Liouville 
integrability of the classical system, i.e. the existence of n independent integrals in in- 
volution, was already proved by Moser |IVIc| using a Lax pair formulation. For a partial 
gener alization of this result to the other classical (i.e. non-exceptional) root systems, see 
| OPl| ] and [@. 

Just as the hypergeometric PDO, Macdonald's difference operator for R = An-\ is 
related to certain known quantum systems of n particles. The systems of interest were 



originally introduced as a relativisti c ge neralization of the Calogero-Moser systems |RS| 
(classical) and [Rl] (quantum). (See |R2 for a survey and connections with certain soliton 
PDE's and exactly solvable quantum field theories). The relativistic generalization of 
the Calogero-Moser system is (quantum) integrable too: explicit formulas representing n 
independent commuting integrals are presented in |RS|, R1|. At the quantum level these 



integrals are AAO's, which are related to Macdonald's yl„_i-type difference operators E,^^ 
via a similarity transformation |K1|. The operators E^^^ (associated with the fundamental 
weights uor) generate an algebra of commuting AAO's which have Macdonald's j4„_i-type 
polynomials as joint eigenfunctions. 

Also for root systems other than An-i one expects that there exist algebras of commut- 
ing difference operators that are simultaneously diagonalized by the Macdonald polyno- 
mials. The purpose of the present paper is to introduce such difference operators for the 
Koornwinder polynomials. We will show that via limit transitions and/or specialization 
of parameters, this also leads to the corresponding AAO's for those families of Macdonald 
polynomials that are connected with (admissible pairs of) the classical series: An-\, Bn, 
Cn, Dn and BC„; the Koornwinder polynomials then reduce to the latter Macdonald poly- 
nomials. By sending the step size of the differences to zero (this corresponds to the limit 
g — > 1) our AAO's go over in PDO's. Thus, we recover the commuting hypergeometric 
PDO's associated with the classical root systems as a limit case. 

Our difference operators constitute a new integrable quantum system of n particles in 
dimension one. In this paper, however, we will not pay much attention to this interpretation 
of the AAO's; instead we will emphasize the connection with orthogonal polynomials. In a 
forthcoming paper the author intends to return to the question of integrability of these and 
related n-particle systems, at the level of both quantum and classical mechanics Q. In 
particular, possible generalization to integrable systems consisting of commuting difference 
operators with elliptic functions as coefficients will be discussed in Q. 

Before outlining the contents of this paper in more detail, let us mention two more 
connections of interest. For special values of the parameters (namely those corresponding 
to root multiplicities), the system of hypergeometric PDO's coincides with the radial re- 
ducti o n of the algebra of invariant differential operators on certain symmetric spaces G/K 
0P2, tIO|. It seems natural to ask oneself the question whether, for special values of the 
parameters, our system of difference operators can be seen in some way as radial reduction 
of certain AAO's connected with quantum homogeneous spaces. Recent results on the 
quantum group interpretation of Macdonald's A^-i-type polynomials [K3, ^ indeed seem 
to point in this direction. However, apart from this special case no relations of this kind 
are known to the author. 

Recently, Cherednik introduced commuting difference operators connected with Knizh- 



nik-Z amolodchikov-type difference equatio 
Ch2|. It is claimed by the author that fc 
donald AAO's. It would be interesting to ii 
of the present paper. Specifically, one woulc 
that are associated with the classical root sj 

The paper is organized as follows: in Sec 
r — 1, . . . , n in the (real) variables xi, . . . ,x 
itly) depend, apart from two scale factors, ( 
Di, coincides (up to an irrelevant multiplit 
operator D^^ ; for r > 2 our operators are ni 
reduces to the rth elementary symmetric f 
(Here Di{xj) denotes the one- variable vers 

Section ^ contains the main results of th 
that: I. Dr leaves invariant certain finite-dii 
ii. Dr is symmetric with respect to the L 
being the orthogonality measure of Koornwi 
facts implies that Dr is diagonalized by th 
mention that this method to diagonalize th 
the approach that was originally used by 
eigenfunctions of the (trigonometric) Calog 
the operators is computed explicitly. As a re 
phism between the abelian algebra general 
in n variables. For Di the discussion in thi 



already obtained by Koornwinder | K4 1 . Thf 
stricted to the case r — 1) and that of | K4 
and some asymptotics) we avoid certain (ra 
leading to the triangularity and the eigenvs 

In Section ^ we discuss the transition tc 
ometric PDO's related to the root system 
q ^ 1, the eigenfunctions converge to the 1 

In Section ^ we study various special casei 
introduce a limit transition leading to the / 
from a physical viewpoint. Specifically, we s 
AAO's can be obtained from the Koornwi 
respectively. This novel transition can be aj 
it enables one to view the A„-i Jacobi pol; 

In 5.3 and 5.4 we show how the Macdon 
root systems can be obtained from the Kooi 
of the parameters. In contrast to our ac' 
contained, this involves various concepts f 
account more accessible by collecting som 
skipped at first reading and referred back t 
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2 Introducing the Difference Operators 

In this section the operators Dr, r = l,...,n, are introduced and their combinatorial 
structure is discussed. 

2.1 The Operator Dr 

In order to write down our difference operators we first introduce some notation. Let Va{z) 
and vi,{z) be the following trigonometric functions: 



Va{z) 



sin a(/i + z) 



(2.1) 



_ sinQ(/io + z) cosa(/ii + z) sina(^o + 7 + cosQ(^i + 7 + z) 

'>^b\z) — : , ^_ ^ ^r; /_ , 3 1 \^-^) 



sin a(7 + 



cos a(7 + z) 



with a, 7 and ^, fis, f^'s, {S ~ 0, 1) complex parameters. For later purposes, it is convenient 
to parametrize 7 according to: 

7 = if5/2. (2.3) 
We form the following multivariable functions using Va and vi, as elementary constituents: 



14. 



j<3' 

W '"a{ejXj + a::ft)t;„(eja;j - Xk), (2.4) 



keK 



with 



J,K C {l,...,n}, Jn/C = 0; £j e{+l,-l}. (2.5) 

The variables Xi , . . . , are assumed to be real. The function V^j-k depends on the index 
sets J, K and on a collection of prescribed signs Sj, j £ J; it serves as a building block 
from which the coefficients of Dr are constructed. The AAO's read explicitly 



Dr 



E 



E 



(-ir 



.7C{l.....ii}, |J|=r 9C,;^C...C,/,=j 



(2.6) 



r = 1, 



with Jo = and 



E^ 



l_d_ 

i dxi 



(2.7) 



Remarks i. \J\ denotes the card inali ty of J and J'^ is the complement set {1, . . . , n} \ J. 
ii. The first summation in Eq. (2.6) is over all index sets J C {1, . . . ,n} with cardinality 
r and over all flippings of the signs Sj £ {+1, —1}, j £ J', the second summation is over all 
strictly increasing sequences of subsets in J: 



§ Ji g J2 g ■ • ■ g Js-l g = J, 1 < s < I J|. 



(2.8) 



in. The exponential exp{—f39j) acts on a f 
variables (complex) shift: 



(e-'"'V)(^i 



Hence, Dr is indeed an anal yti c di fference 1 

IV. The 'hats' in Eqs. (|!6), ( |2^ ) are use 
rather than ordinary (complex) functions o 

V. In the simplest case, i.e. for r = 1, Eq. ( 



Di ^ ^ iii,(£Xj) JJ^ t;a(ea;j 



l<j<Tl 

£ = ±1 



w. It is clear that the operator Dr is i: 
xi, . . . ,Xn. Furthermore, the AAO's are £ 
period: a simultaneous shift of the variable; 
parameters: 

Xj ~*Xj +7r/(2a), j = 1, . 



(see Eq. (|2J) and Eqs. (f^), (|2J)). 

For some purposes it is more convenient t 
such expressions read 



Dr = 



E 



E (- 



JC{l,...,n}. I J|=r 0CJoC;---C;js = J 

ej=±l,jeJ 0<s<r 



(with J_ 



and 



E E 



0<3<r JC{l,...,n}, |. 

e3=±l,i6 



with 



Wl,: 



Wi,o 



E (- 



E 



l<g<p 

1. 



\Iq\=P 



In Eq. ( 2.12 ), D,. is written in terms o f the 
[exp(-/36l)-l]. Accordingly, Jo in Eq. ( pT2 
the fact that the coefficients of the transla 
which does not commute with the translate 

Note. From a physical point of view, one 
particle quantum system in dimension one. 
{S = 0, 1) and Va{ixj lizxt) in the coefficieni 
sible for the interaction between the partii 
field. 
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2.2 Combinatorial Structure and Parameters 

The increment sets Ji, J2 \ Ji, ■ ■ ■ , Ja \ Js-i of the increasing sequence ( |2.8| ) form the 
blocks of a partition of J; the second summation in (2.6) amounts to a sum over all ordered 
blocks. By breaking up Vej-K (Eq. (2-4)) into three parts 



with 



W'"b{ejXj), 

jeJ 

Y\_ ^"fes^j + + Sj'Xj' + 27), 



j<j' 



(2.15) 

(2.16) 
(2.17) 

(2.18) 



one can rewrite Dr (Eq. (|2.q)) as 

JC{l,...,»i}, |.7|=r 



(2.19) 



= - 1 



l<a<r 



Eqs. (2.6), (2.12) and (2.1J) are more compact than (2.19), but the latter has the virtue 
that differ ent parts of the coefficient can be controlled independently. The index set J in 
Eq. (2.1£) will be referred to as the cell. The first block Ji determines the translator; this 



part of the cell will be called the nucleus. Notice that: i. V^jV^j. jc depends on the cell 



J but not on its subdivision in blocks; ii. the product over V^f^j ^-^j 
partition of J, but not on the order of the blocks; Hi. the product over V!^ 



depends on the 

e(J,,\J,,_l);.7\J^ 



depends both on the blocks and on their order. 

The parameters a and (3 are scale factors; a determines the period of the trigonometric 
functions and /3 the complex shift of the translation operators exp(±/36j). Both parameters 
will be taken positive. The parameters fi, fis and fis' determine the relative 'weight' of Va 
and ^6 in the coefficients of the A AO. For instance, Va = 1 for /x = 0; therefore, Va may 
be omitted for /x = 0. The parameters fi, fis and Hs' will be assumed to be non-negative 
imaginary: 

H = i/3g, ^s = iPg6, fi'g = iPg's, (5 = 0,1) 
a,p>0; g,gs,g's>0. 



(2.20) 



Notice that the above restrictions on the parameters guarantee that: i. exp(±/36j) yields 
a pur el y im aginary shift; ii. the commuting part of the coefficient, viz. Wj':,r-s (cf. Eqs. 
(2. 13), (2. 14)), is real because Vc{z) — Vc{—z) (c = a,b) for z real. 



If one picks a = 1/2, then Di (Eq. (^ 
live constant, with Koornwinder's differen 



parameters used in [K4| are related to ours 



-fSgo 



Note. The parameters g, gs and g'g can b 
constants that determine the strengths of t 
setting g — (i.e. ^ — 0) yields a syste 
external field. 



2.3 g = 0: Reduction to Ran 

By setting g — 0, t he co mbinatorial strud 
coefficients in Eq. (2.1£ ) no longer depenc 







= 1). It will be shown 



elementary symmetric function of the follow 



Diixj) = ^ Vb{exj) le 



(For n — 1, Di (2.1C) coincides with Di{xi 



By summ ation of all terms in Dr which c 
Ji = /, Eq. ([2. 191) reduces to 



where 



with 



Dr 



Co 



E 

JC{1, .-.,«}. |J|=r 

ej=±l, jeJ 



1, 



1< 



E 

PlH \-Ps=p 

Pj>l 



To verify this, think of Np^s as the numb 
objects can be distributed over s distinct si 

Lemma 2.1 

One has 
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Proof 

The above interpretation of Np^s leads to the recursion relation 



E 



s— l<q<p— 1 



p - q 



N^,s-i, (p>s>l), 



Substituting (2.2^) in (2.25) yields a recursion relation for Cp 



E 



0<i2<p-l 



whose unique solution is (2.27) 



Np,i = 1. 



(2.28) 



(2.29) 



Definition 2.2 Let ti, . . . ,t„ belong to a commutative algebra. The rth elementary sym- 
metric function Sr of ti, . . . ,tn is defined as 



Sr{ti,...,t„)= ^ n*^' r = l,...,n. 



JC{l,..-,Tl} jsJ 

|J!=r 



(2.30) 



3 Simultaneous Diagon 

In this section it is shown that Koornwinde 
tions of Di , . . . , Dn ■ We prove that the dif 
eigenvalues. As a result, we obtain an expli 
the abelian algebra generated by Di, . . . ,j 
For convenience, we will put a = 1/2 from 



3.1 Trigonometric Polynomi 

Let A = C[exp(i2;i), . . . , exp(i2;„)] be the 
n-dimensional torus 

T =]I 

A is spanned by the Fourier basis {e"^} wit 



\x) 



Let W be the (Weyl) group of permutati 
(so W S„ K (Z2)"). The subalgebra 
polynomials on T is spanned by the basis \ 



Proposition 2.3 If g = 0, then 

Dr = SriDiixi),...,bi{x„)), r = l,...,n 



(with Di{xj) defined by Eq. {2.2i )) 



Proof 

Substituting (|2.27]) in (Eii) yields 



JC{l,...,Ti}, |J| = r ICJ 
Sj=±l,jeJ 0<|/|<r 



Using Eq. (2.16), one completes the proof of the proposition: 



Dr = 



JC{l,..-,Ti}, |J| = r jSJ 

e,-=±l, jSJ 



E n^i(^^) = 5.(61 (xi),...,£)i(x„)). 



JC{l,..-,n} jSJ 



(2.31) 



(2.32) 

(2.33) 
(2.34) 



Note. Proposition is in accordance with the previously noted fact that for 3 = the 
particles of the quantum system become independent. 



mx 



(.)= ^ 



X'eWX A'SS„A 

with ~ denoting proportionality and dt 
P+ = {A e P I Ai ; 
The lattice V can be partially ordered in 
Definition 3.1 [partial ordering of V] 
(VA,A'eP): A' < A iff 



The above ordering induces a partial or 
each dominant weight A G we associa 
highest weight A: 

AY = span{mA'}(A' 
Occasionally we will also use the notation 



E 



l<j<n 
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3.2 Triangularity 

In this subsection it is shown that Dr maps the highest weight spaces into itself. 
Definition 3.2 A linear operator D : — > is called triangular iff 

D{AY) C AY, VA g P+. (3.8) 



One can rewrite Eq. (3.8) in a more illuminating way: 
(VAGP+): Dmx= } ' [D]x,y my , 



E 

\'eV+, A'<A 



with [D]), 



(3.9) 



(i.e. [-D];^ = if A' ^ A). In order to prove that Dr is triangular, we first need to verify 
that the operator maps A^ into itself. 

Proposition 3.3 (invariance of A^) 



DriA'^)CA'^, 



(3.10) 



Proof 

Acting with Dr Eq. (2.1t) on a monomial m\ (3. J) yields the following W-invariant 
trigonometric function on the torus T: 



{Drmx){x) = 



E 

JC{l,.-.,n}, |.7|=r aCjjC...C = j 

Ej=±l, jeJ l<s<r 



(3.11) 



tl<s'<s J 



with 



eeJ 



E' 



(3.12) 



({ei,...,e„} denotes the standard basis of K"). The r.h.s. of (3.11) is rational in the 
exponentials exp( ixj), j = 1, . . . ,n. In order to prove the proposition, we need to show 
that Drm\ (3.11) is actually a polynomial in exp(ia;i), . . . , exp(ia;„). Since the r.h.s. of 
(3.11) is symmetric in x^, . . . ,Xn it suffices to verify that Drm\, viewed as a function of 



xi, is free of poles 

As a function of x\, the terms in ( |3.1l| ) may have poles caused by zeros in the denom- 
inators of the coefficients of the AAO. These poles are located at (cf. Eqs. (2.16)-( p.l' 
and (|2T|), ^^): 

a;i = mod vr, 

= ±7 mod TT, 

i = 2 n 

2 



(3.13) 



• 1 '"1 

. ,n. 



= isj ± 27 mod 27r, j 

From now on the parameters 7, ^, /if, ^'g (5 — 0, 1) and the remaining variables X2, ■ ■ ■ ,x„ 
are fixed in general position. Specifi cally, we choose these parameters and variables such 
that the poles in the terms of (3.11) are simple. 



The residue at xi = vanishes becau; 



Xj vanishes because (3.11) is invariant unc 



Furthermore, because Dm\(x) is even in Xj 
over half the period (cf. Remark vi, Sectioi 

Xj — > Xj + TT, J = 1, . . . , n <^ 

(with |A| = Aj), we need only show t' 



27 



type I: xi = -7. 



The only terms in the r.h.s. of (3.11) that • 
corresponding to cells J with 1 £ J and ei = 
a cell J and choose a configuration of signs 
the total residue at xi = —7 in the sum of al 
J, with the signs prescribed, is zero. One n 
this is because the general situation can be 
flipping of signs of t he va riables Xj, j £ J. 
the sum of terms in (3.11) corresponding tt 
regular at xi = —7. 

type II: xi = — Xj — 27. 
The proof of this case is very similar to th 
cont ribute to t he residue at x i = — Xj — 27 
(|t|) and ( ^.17| )). Lemma [a!^ of Appendix . 
correspond to a flxed cell J with the signs e 
Again the general case (corresponding to ai 
be obtained by an appropriate flipping of t 

We conclude that Drm\ is a W^-invari; 
poles. Consequently, DrVfix must be a pol 
the proposition. 



Proposition 3.3 says that DrUix is a W-r 



must be a finite linear combination of mom 



(VA G P+) 



Drtnx = 



A'6 



with 



V + r = {A' G 7 

In order for Dr to be triangular one must 1 
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(cf. Eq. (3.9)). We shall prove this property by studying the asymp totic s of { D^ m x)(x) 
for Im Xj — > —00. The following limits will be useful (cf. Eqs. (2.1), (2.2) and (2.20)): 



lim Va{z + ieR) 
lim vt{z + ieR) 

R — voo 

(with e = ±1). 

Proposition 3.4 (triangularity) 



e/39/2 



e/3(90+91+9o+9l)/2 



(VAeP+): Z).(^r)C^r, r = l,...,n. 



(3.19) 
(3.20) 

(3.21) 



Proof 

Fix an r € {1, . . . , n} and A G "P"*". Let u)h = X]i<j<j. Cj (the fcth fundamental weight) 
and introduce (cf. (5.17)) 

Mx,r;k=Yaaji{{\' ,LOk) \ \'e'Pt^r}- 



(3.22) 



To derive a contradiction, assume ( 3.18 ) does not hold; i.e. assume that there exists a 
k £ {1, . . . ,n} such that 



Mx.r:k > 



i\,u;,) U J2 

\ i<j<k ) 



(3.23) 



Now it is easy to verify the asymptotics 

mA/(a;-ji?ajfc) ~e^'^''"''=^ ^ e"" {x), R ^ +oo, (3.24) 



with 



Wy.k = {weW\ {w\',Uk) = {\',Uk) }, 



(3.25) 



so using ( p.l6| ) we obtain 

' {Drm\){x - iRcjk) = 



lim e 



E 



(t) 



,A"ew,,^,(A') 



(3.26) 

On the other hand, Eq. (B^) combined with the limits (3.19) and (3.20) entails the 
following asymptotics for (3.11): 

{Drmx){x~iRLUk)^0{e^'-^-'^'''>), R ^ oo. 



Consequently, 



lim e 

R^oa 



(3.27) 
(3.28) 



(because of inequality (3.23)) 



(a;) (3.2) corresponding to di fferen t wei ghts A " G V are linearly independent. Hence, 



The matrix elements [-Dr]A,A' in (3-26) are non-zero (by definition), and the exponentials 



by comparing the r.h.s. of Eqs. ( p.26| ) and (3.28) one arrives at the desired contradiction. 

□ 



3.3 The Spectrum 



By extending < (Definition 3.1) to a linear 
ordering), it is easy to see that the triangu' 



Drmx 



E 



has as consequence that the elements on th 
values of the operator £>,■ : — *■ . (It 
these eigenvalues are semisimple). The pu 

[■Dr]A,A. 

Let y G ]R" be a fixed vector subject to t 

yi > y2> 



By combining the asymptotics (cf. (3.24)) 



mx'iiRy) 



with Eq. (3.29), one finds (use A' < A ^ J 



[Dr]x.\ = lim e 

R — >oo 

In the next proposition, we will evaluate tl 

Proposition 3.5 (eigenvalues) 

One has 



with 



[Dr]x,x = 2'' £;,,„(ch/3(Ai + pi), . . . , 



Er^n 1 ; ■ ■ • ; in i pr ; • ■ • ? Pn ) 

/ 



E (-1)'- 



E n 

jc{i,..-,n} je. 

L7|=s 



and 



p] = in- j)g + {go + gi + 



Proof 

Using Eq. (|2.13|) we obtain 



{Drmx)(x)= J2 E 



0<S<I- ./C{l,.-.,n}. |./| = ^ 



7 



with Vs.j: K and W^/,p defined by Eqs. (2.4) and (2.14), respectively. In order to compute 
the limit (3.32), w e firs t derive some preliminary asymptotics: 



I. mx (cf. Eq. (3.24)) 



mx{x + if5eej)\:c=iRy ~ 6 



i? — > 00. 



II. Ve. 



From (3.19), (3.20) one deduces 



fim V^j;j'i\x=iRy = e' 

J? — »CX3 



/3[9iV..7 + (90 + 91+9o + 9l)A^ejl/2 



with 



|{j e J I = +1}] - |{j G J 1 = -1}| 



(3.37) 



(3.38) 



(3.39) 



and 



2\{j,keJ\j<k, =+l}| -2\{j,keJ\j <k, = -l}| + 

2 \{j eJ, ker\j<k, e, = +i}| - 2 |{j eJ, ker\j<k, e, = -i}| 

2 |{j G J, fc = 1, . . . , n I j < fc, = +1}| - 
2\{j eJ, k^l,...,n\j <k, Ej ^ -l}\ 

2^(n-j>,. (3.40) 



Consequently, 



lim K 



(3.41) 



with pj defined by (3.35) 
in. Wi,n 



Using (2.14) and (3.41) it is not hard to see that fim_R_oo Wi ,p\x=iRy exists and depends 
only on the cardinality of / and on p (but not on the number of variables n). We define 



Notice that (cf. Eq. (2.14)) 



F|7|_p = 2 lim Wi,p\x=iRy. 



Fm.o = 1, 



(3.42) 



(3.43) 



After these preliminaries, we are now ready to compute limit (3.35). Substituting ( p. 3' 
in (3.32), and making use of (3.3'i), (3.41) and (3.42), we obtain 



[Dr]x.x = 2'- ( n +pAf« 

.7C{l,...,n}, lJl=s \j<^J J 



(3.44) 



It remains to c alcu late Fm^p, 1 < p < m < r 
zero (see Eq. (2.6)); one obtains, therefore 
A = 0: 



j2 n^^^^-"^ 



JC{l,...,Tl}, |J|= = 

0<s<r 



For a fixed number of variables n, this y 
same coefficients Fm,p occur, now with '_ 
for n' — 1, . . . ,n, and making use of Eq. ( 
equations in the n(n + l)/2 variables Fm,p, 
In order to solve this system it is convf 
Notice that pj (unlike pj) does not depei 
substitution the n{n + 1) /2 equations beco 



JC{1 n'}, |.7| = s XjeJ / 



with condition (3.43). In Lemma 



B.l 



of A 



has a unique solution: 



Fm,p = (-1)'' 



Subst ituting (3.47) in (3.44) and using pj 

(Esi). 



For g = all components of the vector p 

pi,P2,...,Pn = (50 

Then, one can rewrite the above expressic 
symmetric functions (see the remark follow 

[br]x.X = 2'- Y li^'^ 

jc{i,.-.,n} jeJ 

\J\=r 

= 2" Sr (ch/3(Ai + p] 
This equation is in agreement with Proposi 

3.4 Symmetry 



In Ref. [ |K4[ the following weight function c 

A{X) = Y[ da{Xj + 

l<j<k<n 
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with 



(•g-/3ffgiz.g~,3)^ 



d+iz) = 



c = a,b 



(e ;e 



The so-called q-shifted factorials are defined in the usual way: 



aq 



(3.5f) 
(3.52) 

(3.53) 



(3.54) 



Notice that the conditions on our parameters, viz. ( |2.2C| ), guarantee that the infinite 
products in E ns. ( js^S^ ) and (3.5S) converge. Recall that in order to compare our formulas 
with those of [K4|, one has to reparametrize according to Eqs. (2.21) and (2.22). 

Let L'^{T , Adx) be the space of W- invariant functions on T that are square integrable 
with respect to the measure Ada;. We define 



if, 9) A 



fgAdx, f,geLir{T,Adx). 



(3.55) 



The space of W-invariant polynomials is a dense subspace of L'^{T , Adx). The 
purpose of the present section is to show that the A AO's Di, . . . , Dn are symmetric with 
respect to (•, •)a. We need the following lemma: 

Lemma 3.6 Let a = 1/2 and z G M; furthermore, let the parameters be subject to condi- 



tion (2.2L ). Then the functions di^"* (c = a,b) satisfy the following first order difference 
equations: 



d+: 



ii. dc : 



Proof 



dt{z + ip) = e-''''^va{z)dt{z), 

d+{z + i(3) = e-^'''"'+'>^+'>'o+^'^^/^ Vb{z)dt{z)- 



Vc(z — j/3) dc{z — if)) — Vc{z) dc(z), c = a,b. 



Eq. (|3.56|) is an immediate consequence of definition (3.5^ 



dt{z + if3)/dt{z) 



(1 - e-^'^e'") 
(1 - e'^) 



-0a/2 



Vaiz). 



(3.56) 
(3.57) 

(3.58) 



(3.59) 



Eq. (3.57) can be reduced to the former case by observing that d'^{z) factorizes: 



dtiz) 



(e";e-'3). 



(- 



(e-/3soe»2;e-'')oo (-e-'^fie'^; e-'')o 



(- 



(g-/3K + l/2)giz. g-/3)^ (_e-/3K + l/2)g>z. g-/3)^ 



(3.60) 



ii. Using Eqs. (3.56) or (3.57), respectively 



Vc{z — if3) dc{z ~ if3) 



Vc 

dt 



Part ii. of the above lemma leads to the fo 
Corollary 3.7 One has 

We now arrive at the main result of this 



proof hinges on relation (3.62). 
Proposition 3.8 (symmetry) 

{Drmx,mx')A = {mx,L 

Proof 

First consider the following contour Integra 



W. 



J<',r-s VsJ-J' 



with Vej-K and as in ( |2.4D and ( p.l4| ) 

place over the closed contour 

Cj — [— TT, tt] U [tt, tt — iSjP] U [tt — ie 

Let all parameters and the variables Xk, i 
the integrand has simple poles inside Cj d 
A(a;). However, one easily verifies that ai 
a zero in A{x); similarly, pol es in side Cj 
Vej-ji. Consequently, inte gral (3.64) vanish 
the contributions to (3.64) which are due i 
respectively, cancel each other because the i 
upshot is that, when integrating the integr; 
integration path to [— tt — iejl3, tv — iejjS] w 
Armed with this conclusion and Eq. ( ^ 
can be written 



{DrTUx, mx')l 



Eq.(2.13) 



E 

|J| = s,0<s<r 
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Deformation of the integration paths of Xj, j £ J, from [— tt, tt] to [— tt — iejl3,-K — iSj/S], 
followed by a change of variables Xj Xj — i/Sej, j £ J, yields 



{Drmx,mx')A= ^ 



\J\=S,0<3<T 



Wj'^,r-a mx (e*^*^-^ mx' V^j-j-^ A)dx. 



(3.67) 



Using Corollary 3.7 and the fact that Wji^r-s is real (for parameters subject to ( |2.2C| )) 
entails 



{Drmx,mx')A = ^ / mx Wjc^r-s Vej;j<' {e'l^^^J nix') Adx 
|j|=.,o< = <.-^T 

= {mx, Drmx>)A- 



3.5 Diagonalization and Commutativity 



(3.68) 
□ 



If g — go — gi — g'o — g'l = 0, then da — dt — 1 (see Eqs. ( |3.5]| ),(3.52) and (3.6C)), and thus 
{■, •)a reduces to the inner product on T with respect to Lebesgue measure (A = 1). The 
basis of monomials {mx}xev+ orthogonal basis of L'^{T , dx). For arbitrary parame- 
ters however, the orthogonality of the monomials with respect to (■, ■) a no longer holds. By 
subtracting from mx the orthogonal projection of mx onto spa.n{mx'} (^x' ev+ , \' <\) 
one obtains an alternative basis {px}xev+ of ■ This is the basis of Koornwinder poly- 
nomials. 



Definition 3.9 [Koornwinder polynomials] 

Koornwinder' s polynomial px £ A!x is defined by the conditions 



V\ = mx + 



E 



Ca,a' mx' 



and 



9A,mv)A = 0, VA' £ P+, A' < A. 



(3.69) 



(3.70) 



We now prove that {px}xev+ is a basis of jo int e igenfunctions of Di, . . . , D„. For conve- 
nience, the notation for the eigenvalues Eq. ( 3.33 ) is sometimes abbreviated by putting 



£;.,„(ch/36'i, . . . ,ch/30„;ch/3p„ . . . ,ch/3p„) — ^ Er^^O), (6 £ M"). 



Theorem 3.10 (eigenfunctions) 

Dr Px = Er.n{\ + p)px, 

(with r = 1, . . . ,n and p = (pi, . . . , pn), see Eq. (3.31 )). 



VA £ 



(3.71) 



(3.72) 



Proof 

It follows from (3.69) and Proposition 3.4 t 



{Dr - [Dr] 

On the other hand, one has (use the propoi 



{{Dr ~ [Dr]x.x)px,mx')A (p^ , {D^ 



Combining (3.73) and (3.74) entails Dr px 
to complete the proof. 



In Appendix C it is shown that if a diffe 
then all its coefficients must be zero. Comb 
commutativity of the A AO's: 

Theorem 3.11 (commutativity) 

The operators D\, . . . , Dn mutually commu 



Proof 

The polynomials {px }x9V^ 
D\,...,D„ (Theorem 3.10). Hence, it is 



form a basis 



on ^ . In other words, the commutator 
vanishes on A!^ : 

[br,br>\{A'^')=Q, 

It now follows from Proposition |C.l| in Appc 
are identically zero. 



Consider the real algebra of difference of 

D = M[i 

It is clear that D is an abelian algebra (T 
simultaneously diagonalized by the Koorn\^ 

Theorem 3.12 (D 5^ Mfch/S^i, . . . , ch/36'„]' 
For each symmetric function S{9) £ IR[ch/36 
operator £) £ D such that 

Dpx = S{X + p) 

Proof 

Er,n{9) (3.71) is a linear combination of ek 

Er,n{e) = Sr{ch(39l,...,chp6 
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l.d. stands for terms of lower degree in chf36j, j 



The elementary symmetric 



functions form a set of algebraically independent generators of the sym metr ic algebra 
(this fact is the 'fundamental theorem on symmetric functions', see e.g. |M1|). Hence, 
Eq. (3.78) implies that the same is true for the functions Er,n{0): every element in 
R[ch/36'i, . . . ,ch/36'„]^" ca n be written uniquely as a polynomial in Er,n{6), r — 1, . . . ,n. 

Now we use Theorem to conclude that for every symmetric function S{9) there 
exists a difference operator D G D such that Eq. (3.77) holds. That such a difference 



operator D is unique follows from Proposition C.l (Appendix C) 



The symmetric functions S{6) G M.[chf36i, . . . , ch/36'„]'^" separate the points of the wedge 

{e eR" \ ei>e2> ■ ■■ >0n> o}. (3.79) 

To see this, first notice that ( — l)"^''5'r(ch/36'i, . . . , ch/3&„) is the coefficient of i/""'' in the 
characteristic polynomial det(T — i/I) of the diagonal matrix T = diag(ch/3^i, . . . , ch/3S„). 
Conse quently, the values of Sr{chl36i, . . . , ch/36'„), r — 1, . . . ,n determine 9 in the wedge 



(3.79) uniquely. This fact combined with Theorem 3.12 can be used to prove the orthogo- 
nality of the basis {p\}: 



VA,A' G P+, A / A'. 



Corollary 3.13 (orthogonality) 

(pa,Pa'>a = 0, 
Proof 

Let A, A' G ; and A 7^ A'. S ince the symmetric functions in IR[ch/36'i , . . . , ch/3&„ 



(3.80) 



such 



separate the points of the wedge (3.79), there exists an S{d) G R[ch/3di, . . . , ch/3^„] 
that 

S(A + p)/S(A'+p). (3.81) 
Hence, by Theorem 3.12 there exists a D G D for which p\ and px' are eigenfunctions 



corresponding to different eigenvalues. But then the polynomials px and px' must be 
orthogonal with respect to {•, ■)a because D is symmetric, cf. Proposition p.q. 



Note. The orthogonality of the basis {px}xev+ already shown by Koornwinder 
His proof exploits the continuity of {px,Px')a in the parameters. 

Corollary 3.14 (self-adjointness) 

Every difference operator £) G D is essentially self-adjoint on C L^{T,Adx). 
Proof 

This an immediate consequence of the fact that every AAO in D acts as a real multipli- 
cation operator on the orthogonal basis {pa}agp+ '^^ ^^(T, Ada;). 



Remarks i. Theorem |3.12| states that the ai 

HC 



Dr 



Er,n(6 



induce a Harish-Chandra-type algebra ison 
ii. Recall that px is defined as mx 1 
spaxi{mx'}\' ^x- Using the orthogonality c 
following recursion relation for p\: 



Px = mx 



3.1 



is extende 



Hi. If the partial ordering 
is possible to orthogonalize the basis [mx] 
Corollary 3.13, the result does not depend 



the ordering: the resulting orthogonal basii 
restrictive property of the measure A dx. 
iv. According to Theorem 3.11, 



the Hi 

integrable n-particle system (cf. the note e: 

4 (3^0: The Transitu 
metric PDO's 

By sending the step size j3 of the differences 
pergeometric PDO's associated with the rot 
{px} converge to the BC„-type Jacobi polj 
In this section we will make the dependent 
all objects of interest, e.g.: -Dr,/3, A/3 and p 

4.1 Eigenfunctions 

Consider the following weight function on 1 



Ao{x) = 



l<j<k<n 



n 

n I 

l<j<n 



I sinQ(a;j + Xk) 



|2go 



and let {•, •)ao be the inner product on L 
introduces W^-invariant polynomials on T 
Definition 3.9); these are the BCn-type Jac 



Definition 4.1 [BCn-type Jacobi polynon 
The Jacobi polynomial px.o G A^ is definet 

PA,o = rriA + 



11 



and 



(pA,o,mv>Ao = 0, VA'gP+, A'<A. 



(4.3) 



Remark. Usually Ao is written in a slightly diflterent form, which emphasizes the relation 
with the root system BCn- Use 



|sin(a2;j)| ®° |cos(Qa;j)| ~ |sin(a2;j)| " |sin(2aa;j) 



2fei 



(4.4) 



with go = fco + fci and gi — ki, to compare ( [4.1[ ) with the usual expression for the i3C„-type 
weight function. 



We need the following convergence result from |K2| to connect pa.o with Koornwinder's 
polynomial Pa,/3, /3 > 0: 



i'=2-fcl 



Lemma 4.2 Let fci, fc2 £ M and q g]0, 1[. Then 

uniformly for z in compacts of the punctured disc {z (z C \ \z\ < 1, z ^ 1}. 
Proof 



(4.5) 



See Proposition A. 2 of Ref. |K2, Appendix A] 



It is immediate from Eqs. (3.52), (3.6C) and the above lemma that 
lim daA^) = I2sin(2/2)|^», 

;3-.0 



lim db,piz) = l2sin(2/2)|^<™+»o^ \2cos{z/2)f^^^+'''^\ 



(4.6) 
(4.7) 



uniformly for z in compacts oflvr, 7r[\{0}. 

Consequently, for /3 — ^ 0, A^g (3.5C) converges to a weight function that is proportional to 
Ao: 



lim Art = c Ao, 



with 



go = go + g'o (> 0), 

and c = 4"("-i'f 4"(»+si) (> 0). 



gi ^ gi+g'i (> o) 



Proposition 4.3 One has 



lim pA,/3 = Px,o, VA e P 



and 



3a.o,Pa'.o)ao =0, VA,A' eP+, A/A'. 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



Proof 

Using (|3.8^), (4.^ ) and an induction argume 



satisfies a rec ursion relation of the type (3.8 
and Eq. (4.8) to conclude that the resultir 
{■, ')ao- But then the recursion relation fo: 
the orthog ona l projection ( with respect to 
Definition 4.1 one obtains (4.1C). 



Note. The limit (3 —> corresponds to tl 
polynomials the q 1 limit to the Jaco 



studied in [M2 



was proved in 



M2 


, fc 


in 


HI 



HI I (again for arbitrary roo 



4.2 Eigenvalues 

The purpose of this subsection is to inves 
(Proposition 3.5) for P 0. 



Proposition 4.4 One has (for r — 1, 

-2r jp i^-uaa , _ ^ ch/S^n ; ch/j 

9^/2; p?/2,.... 

Proof 



lim (3'^'' ErJchpei, 

13-0 

= Er,n{9i/2, . . . , 



Our proof of E q. (4.12) hinges on a recursi 
to Lemma B.2 of Appendix B: 



Er,n{chl36i, chfJOn; chPpr, 

(ch/36»i,...,ch/3e„_i 

with the convention 

Eo,n = 1, 

We divide Eq. ( [l.l3| ) by and then use i 
lim /3"^''Sr „(ch/36li, . . . , chpO, 

= (^n/2 — Pn/2) -Br-l,n-l 



+Er,n-i{di/2, . . . ,9: 



2 



Now we use again Lemma 

(Ei3)). 



B.2 to conclude 
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4.3 Operators 

Expansion of -D^./a in P yields a formal power series of the form 

Dr.fi = ^ I)<'"^/3'". (4.16) 

The coefficients are polynomials in the partials 6j, j = l,...,n; this means that 

these coefficients are PDO's. We define t he le ading differential operator Dr,o of Dr,f3 as 
the first nonzero coefficient in expansion (U.ld): 



(4.17) 
(4.18) 



Definition 4.5 [leading PDO] 
Let 

m.r = mm{m G N | 7^ } 

(with Dr^^ defined by expansion ([l.l6| )). Then, 

is called the leading PDO of D,.)}. 

The BCn-type Jacobi polynomials are joint eigenfunctions of -Di,o, 
Theorem 4.6 One has nir = 2r and 

A-,0 PX.O ~ Er,n ((Al + pi) , . . . , (A„ + pn) \pr, - ■ ■ , pn) P\.0, 



VA e 



(4.19) 



uiit/i r — 1, . . . ,n and p as in (3.3t 



Proof 

Consider the eigenvalue equation (3.72): 

Dr,i3 (p\,p) = 2'' £,.,n(ch/3(Ai + pi), . . . , ch/3(A„ + p„);chf3pr, . .., ch/3p„) pA,,3- 



(4.20) 



First, appl y Tay lor's theorem to the l.h.s. of Eq. ( [l.2C| ) and make use of Definition 4.5 
and Limit (4.10) to conclude that 



Dr,l3P\,0 = Dr,OP\,0 P"^'' + o{f3"^'' ) . 



(4.21) 



Next, use Pr oposi tion and Limit ( [4.10| ) to derive the asymptotic behavior for ^ of 
the r.h.s. of (|^): 



2'' Er,n{ch/3{Xi + pi), . . . ,ch/3(A„ + p„); ch/3pr, . . . , ch/?p„) px,i3 = 

Er,n{i\l+Plf,...A^n+Pnf;Pr,---,pl)pX,0P^'' +o{(3^^). (4.22) 



By Definition 4.5 and Proposition CI it is possible to pick a A G P"*" such that Drfl Px.o 7^ 0, 
so rUr > 2r. It is not difficult to see that there also exist A G such that 

Er,n ((Al + pi)^, . . . , (A„ + p„f; pi,..., pi) / 0, 



Corollary 4.7 



Dr 



lir 



Explicit computation of ( [l.23| ) for r = 1 yi( 
Di,o 



l<j<n 



fc <n 

{go cot(aa 



— 2iag y { cot a 

l<j<fc<n 

-2^a J2 

l<j<n 

(with a = 1/2 and go, gi as in Eq. ( |4.9| )). 

An immediate consequence of Theorem | 
Dr,o, r = 1, . . . ,n. 

Corollary 4.8 The differential operators 1 

Let Do = M[i)i,o,...,-Dn,o]- The algebra 
simultaneously diagonalized by the BC„-ty 

Theorem 4.9 (Do = R[el, . . .,elY'^") 
For every symmetric polynomial S{9) G M| 
operator Do G Do such that 

Do Px.o = S{X + p 



Corollary 4.10 (self-adjointness) 

The PDO's Do G Do are essentially self-ad 



The proofs of Theorem 4.9 and Corollary ^ 
and Corollary S.l'^, respectively. 



Remarks i. The operator -Di,o (4.24) c 



so rUr < 2r. This entails m,- = 2r and (4.19) 



hypergeometric PDO associated with the 1 
with the usual notation which emphasizes tl 
ta.n{axj) from (4.24) by means of the rela 
the remark under Definition 4.1). 

a. The existence of an abelian algebra c 
to M.[9'f, . . . , Ol]^" via a Harish-C handra-t yi 
shown by Heckman and Opdam [HI, O^ . 
consider arbitrary root systems. 

Hi. For an arbitrary root system R, th( 
ing generalized Calogero-Sutherland s yster r 
PDO via a similarity transformation | 0P2t 
amounts to the conjugation with A^'^ , wh 
Plancherel measure with weight function A 
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H = A,y'l)i,o A"''" + Eo 



Xk) + sin ^ Q(a;j - Xk)} 

l<j<n l<j<k<n 

+ ^ {go(ffo - 1) sin~^(axj) + - 1) cos~^(aa;j)} (4-26) 



l<j<n 



with Eo = 4q^(p,p) and pj = (n - i)s + {go +gi)/2, cf. (3.3£). 

Corollary i.i can be interpreted as the quantum integrability of the BCn-type Calogero- 
Sutherland system. For arbitrary root systems integrability follows from [ |Hl| , p2| (cf. 
Remark ii). 

5 Special Cases Related to Classical Root Systems 

By limit transitions and/or specialization of the parameters, the operators Di, . . . , D„ re- 
duce to commuting AAO's, which are simultaneously diagonalized by Macdonald's poly- 
nomials. Such difference operators are obtained for all Macdonald families associated with 
(admissible pairs of) the classical root systems: A„_i, Bn, C„, Dn and BCn. 

Note. Most results in this section have an obvious counterpart for /3 = (which amounts 
to g = 1). 

5.1 Preliminaries 

First, we o utlin e very briefly some of the main points of the constructio n pr esented by 
Macdonald |M2j. A more detailed summary of his results can be found in |M4| and [K4]. 
(For our purposes, especially the second summary is useful). Here, we only want to intro- 
duce some terminology which facilitates clarifying the connection between the preceding 
sections and Ref. [M2|. For general information on root systems the reader is referred to 
e.g. 1^ |[. Ahhough most of the remaining part of the paper should be accessible without 
a detailed knowledge of root systems, a glance at the 'planches' in Bourbaki Q might be 
of some help . 

Ref. |M3] uses the concept of admissible pairs of root systems. The pair {R, S) is 
admissible if 7? and S are root systems (assumed irreducible) such that S G R is reduced 
and generates the same Weyl group as R. Let V be the real vector space spanned by R 
and consider the torus T_r = V/{2-kZR^). Let be the dominant cone of the weight 
lattice of R (which equals the character lattice of T_r) and let denote the algebra 
of W-invariant (trigonometric) polynomials on T_r (this algebra is isomorphic to the W- 
invariant part of the group algebra over the weight lattice). To every admissible pair 
{R, S) Macdonald associates a weight function Aj^^g) on T_r and finds a corresponding 
orthogonal basis {p\,{r,s)} xev^ -^'r ■ Furthermore, he introduces difference operators 
Da associated with the so-called (quasi-) minuscule weights cr of S^. These operators are 
diagonalized by the basis of Macdonald polynomials {p\,{R,s)}xev+' 

We will show that additional AAO's for Macdonald's polynomials arise as special cases 
of Di, . . . ,D„ (2.6). This leads to difference operators associated with the fundamental 
weights of S'^, for every admissible pair consisting of classical root systems. These AAO's 



generate an abelian algebra Dj^^g) of differ 
onalized by the basis {^'A,(i^,s)};^gp+ • The 
W-invariant part of the real group algebra 

Remark. If all roots in R have the same lei 
S — R and there exists only one admissible 
then there are several possibilities for the j 
are six such possibilities; these correspond 



5.2 The Root System A„_i 

Let Dr.iead cousist of thosc ter ms in Dr 
exp(-/30,), j = l,...,n (cf. Eq. (El|)): 



Dr,lead ~ ^ ] Vj;. 



JC{l,...,n} 
|J!=r 



One picks up these leading terms via the fc 



with 
Let 



Dr,i„ad = lim e 

R^oo 



\<i<k<n 



Conjugation of Dr lead with A+ results ii 
Kg)): 



D'r ~ A+ _Dr,iead A I 



n n '^"(^^ 

Jc{i,..-,n} \ jej 
\.j\=r \fce./= 



It is clear from Theorem i.ll and Eqs. ( 



commute. For r = n, (5.6) reduces to an of 
coordinates: D'„ = exp {^—f3{9i -!-■■ + ^„)) 
parts: 

D'r = D 

The first part causes a translation x x + i 
(up to a multiplicative constant) with Maci 



Eujr — Cr Dr,An~l : 



Cr 
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The operator E^^^ is assoc iate d with the rth fundamental weig ht ujr of the root system 
An-i. The parameters in [M2| are related to ours via Eq. (2.21). 



Notes i. For r = 1, Eq. (5.6) reduces to 



(5.9) 



l<j<n Xki^i 



li. After transformation to Lebesgue measure, the operator _Dj, goes over in the r- 
th quantum integral Sr [R.1, (Eq. (2.3))] of the relativistic Calogero-Moser system with 
trigonometric coefficients. More precisely, let 



(5.10) 



l<j<fc<7l 



then 



,1/2 



D' A 



-1/2 



1/2 



E 

JC{1 >>} 

J| = r 



a;fe) 



fc6J<= 



/ 

J]^ Ua(Kfc - %) 

\fee J<= 



1/2 



(5.11) 



This relation between the n-particle relativistic CM system introduced by Ruijsenaars 
and Macdonald's difference operators for the root system An-\ was first observed by 
Koornwinder [Kl]. It generalizes the relation between the n-particl e Calo gero-Sutherland 
system and the hype rgeo metric PDO's associated with R — An-i [ |0P2| (cf. Remark lii 
at the end of Section 4.3). 

til. The operators -Dr,A„_i act as a real multiplication on functions that depend only on 
xi + ■ ■ ■ + Xn- The transition D'^ — > -Dr,A„_i can be interpreted physically as restricting 
attention to the motion in the center of mass hyperplane x\ -\- ■ ■ ■ + Xn = Q. 

We show next that the joint eigenfunctions of -D,.,a„_i, i.e. Macdonald's yl„_i-type 
polynomials, can be obtained from Koorn win der's polynomials by a certain limit transition. 
Let nix^i^ad be the sum of terms in rnx ( |3.3[ ) which are of the highest degree in exp{ixj), 
j = l,...,n: 



A'eS„A 



(5.12) 



Recall that according to Definition p\ is a linear combination of monomials of the 
form 

A'£T'+, A'<A 

with Cx \i certain complex coefficients (which depend only on the parameters) such that 
( |3.70| ) holds and ca,a = 1. We set 



PXjlead '■ 



E 



Ca.A' ^y.lead , 



(5.14) 



A'g'P+, A'<A 
|A'N|A| 



Let 

From the asymptotics 



a; = ei 



m\ (x — iRuj) ~ rriA,; 



and Eqs. (5.13) and (5.14), one derives 



P\,leadix) = lim 
R — -oc 



The polynomial px^uad is homogeneous of 
sequently, a translation causes an automoi 
exp(i_R|A|) px,iead- By multiplying px,iead 
ends up with a basis of translation-invariar 



/ \ — i|A| (xiH \-x„)/n 



with 



PX,lead 



-ilX'l 



mx'.A„_i{x) 

As the notation suggests, it will turn out th 
the Macdonald polynomials associated witl 
First we need a lemma. It says that the ! 
(Recall that the eigenvalues of Di are gi-v 
below) . 



Lemma 5.1 Let 



Then, for all A, A' G r+ : 



A > A' 



( with p given by ( 3. 35 ) ) 



Proof 

It is clear from Definition 3.1 that A > A' U 



A = A' + ^ Oj (ej 



l<j<n-l 



with at least one of the Oj 's positive. 



Obviously, it suffices to verify (5.21) fo 
positive. Now for j = n this is immediate, 
convex function: 

ch/3(a; + a) + ch/?( 

if X > y and a > 0. 
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with 



Prop osition 5.2 Let -Dr,A„_i be determined by (5.t )-(5.V) and let p\^A„_i be defined by 
f|j|;. Then 

-C>r,A„_i PA,A„_i = -E'r,A„_i(A + p') (5.24) 

with 



p; = g(7i+l-2j)/2, l<j<n. 



(5.25) 
(5.26) 



(Sr denotes the rth elementary symmetric function (Definition 2.^)) 



Proof 

The operator D'^ (5.(3) is invariant both un der p ermutations of Xj and under translations 
of the form x ^ x + Ruj (with lj as in (5.15)). We u se this and the asymptotics of 



{D'rm\ iead){~iRy) for _R — > oo (with y such that (3.3C) holds) to derive (cf. Proposi- 
tion p^, Proposition 3.4, and their proofs) 



D'r'n^X,lead ~ [D'r]\,y rny,lead 



x'e'P+ , x'<x 
|A'N|A| 



with 



Srie 



(5.27) 



(5.28) 



(The poles at Xj = Xk , j ^ k, cancel because of the permutation symmetry; the condition 
|A'| = |A| in sum the (5.27) stems from the translational invariance of D'r)- 

We will now show that pxjead is an eigenfunction of D'j. (with eigenvalue ( ^.28| )). Con- 
sider the eigenvalue equation (3.72) for r — 1: 



{Drp^){x) = 2( [ch/3(A. 

\l<j<n 



+ p,) - ch/3p,] ] pxix), X£V+, (5.29) 



with Di given by (2.1C). Substitute 



X ^ X — iRjjj 



(5.30) 



and divide both sides of the equation by exp(_R|A|); Sending R oo entails (use 
(El9|),(E2^ and (|5T7|)) 



[c± b'l + C-^^ {b'„) ^b'„-l —C2)p\,lead = 

2 ^ [ch/3(Aj + pj) -ch/3pj] px.iead, 

\1<J<" / 



(5.31) 



^ g-/39("-l)/2g-/3(90 + 5 



C2 



Cl 



l<j<n k^j 



= 2 ^ ch/3p, 



l<j<n 

(To verify equality *, first check that bot 
Xj; consequently, these parts are constants 
the value of these constants by putting x - 
sending R —> oo). 

It is clear from the commutativity of D'l 
on the l.h.s. of (5.31). Therefore, D',.px,i 
corresponding to the same eigenvalue as % 

b'rPx,lead lieS iu Span{ Px> ,lead I A' < A, |A' 

t hat p\, iead must be an eigenfunction of L 
(E2^), 



One obtains the expressions (5.24)-( 5.26) 



0. 



Notice that my^A^^-i ~ 't^a,a„_i iff A' — A G 
can be relabeled by the projection of 
This projection of "P"*" coincides with the c 

The polynomials pa,a„_i can also be rela 

A' - A G Z(ei H h e„). (This follows from 

(5.24) for r = 1, . . . , n — 1). Since the oper£ 
difference operators up to a constant, we ei 



Corollary 5.3 The function px^A„_i coinc 
corresponding to the weight vector X — |A|(f 



as in (2.21)) 



Notes i. The transition pxjead —* Pa,a„_i 
center of mass. 

M. For ^ (i.e. q = exp{-/3) -> 1), tl 



polyn o mials associated with An-i [M2|. It 
(5.14), (5.18) relates the Jacobi polynomial 



-i\X\{xi-\ hXn)/n 



lim e 



lii. Recently, a completely different limit 
BCn to those associated with An-i has b 
with for 5 7^ 1. 
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5.3 The Root Systems S„, C„ and BCn- 

In order to compare our results with [M2| , it is convenient to carry out a reparametrization: 



Ml 



1^1 + V2, 
V2, 



I I 
Ml V2 



(5.36) 



5 = 1,2. 



with (cf. ( |2.20D ) 

= if3ks, ly's = ifJk'g, kg, k's > 
With these new parameters we rewrite vt{z) ( |2.2| ) and d'^{z) (3.5J) (recall also ( [3.6i 

sin a(i^i +1^2 + 2) sin 20(1^2 + z) 



(5.37) 



sina(i'2 + 2) sin(2a2;) 
sin aiu'i + 1/2 + ^ + z) sin 2a{i'2 +7 + 2:) 



and 



4W 



sina(i^2 + 7 + 2) 
(e«(-^2+-).g-/3)^ 



sin 2q:(7 + z) 



(5.38) 



(e*(''2+T-+-);e-'3) 



(e 



2i{7+2). 



-2/3 



(e'(-;+'-2+T+-); e-'3)oc (e2*(''2+T+-); e-2/3). 



(5.39) 



For the following parameters A{x) (3.50) reduces to Macdonald's weight function A(ij_s) 
with R = B„, C„ or BC„ and 5 = 5 or Cfi '• 





Bn Cn BCn 


Bn 
Cn 


' 1 n 
= V2 = V2 = yi 


1/-1 = U2 = 1/1 = 


1/1=1/2=0 


u'l = 

1/2=1/2=0 


l/[=Vl=0 
1/2 = V2 


l/[ = 1/1 
v'2 = V2 



Table. Special cases associated with admissible pairs {R,S). 



The relation with the parameters employed in Ref. [VI2] reads 

, if.L I ih'i 

and 

<? = 



i±2e,- = e 



if S = Bn or S = R = C„ 



211^2 



g-/3/2 S = Cn and i? = B(C). 



(5.40) 
(5.41) 



(In order to verify that for the above parameters, A ( 3.5C ) indeed coincides with the weight 
functions introduc ed b y Macdonald, it may be helpful to compare our expressions with 
Eqs. (3.1)-(3.5) of [K4|, since the latter are rather explicit). 

Next, we consider the Macdonald polynomials associated with A(j{_g). We distinguish 
two cases: 



I. R= {B)Cn 

In this case the torus Th(= M"/(27rZi?'')) o 
invariant polynomials on T_r coincides with 
read 



i^k = ei + ■ 



+ 6; 



(our convention regarding the choice of tl 
of dominant weights V^, which consists of 
k = 1, . . . ,n, coincides with the cone (^ 
2 and 3 of the above table, {pa}a6T'+ redu 
Combined with the structure of the expansic 
with the Macdonald basis {P{r,s),x} x^-p+ ■ 

II. R = Bn 

This case is a bit more complicated because 
the Macdonald torus Tb„(= W /{2-n'lBl) 
weights ujk of Bn are the same as in (5.42), 



l^n = {ei - 

The cone of dominant weight vectors can b 
^i3„ = { A + 

and the algebra of W-invariant polynomial 
mials: 

^2^^ = span{i 



For 5 = (1) the function m\+Stj^{x) (3.3) 

Xj Xj + 2tT =^ f71-A+ 



Combining ( p.46D with the fact that the fui 
zero on the hyperplanes Xj = tt (mod 2n). 

m„„ (x) = 2" 

Thus, we have the following decomposition 



= A 



(with as before). This decomposition i 
on L'^{T B„, ^B^dx) because Ab„{x) is pei 
The situation is now as follows: just a 

sis {pa',{s„,s)}a'6-p+ °^ -^sl splits in pe 

m^^A^ , respectively. By specializing to 
Bn-type Macdonald polynomials that are : 
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B„-type polynomials in m^^A^ can also be expressed in terms of Koornwinder's polyno- 
mials. To s ee this, notice that m^^Ajg^ 5) coincides with a weight function of the type 
(3.50)-(3.53) with parameters 



0, 



S = Bn 
S = Cn 



Ml 



0. 



(5.49) 



By multiplying ma,„ and the Koornwinder polynomials with parameters as in (5.49), we 



obtain an orthogonal basis of rriuij^A ; the latter polynomials coincide with the anti- 
periodic Macdonald polynomials. To be more explicit, we have: 



Pa+ij„,{b„,s) 



Pa 



fio = 1^1, Mo = Ml = m'i = 0, S = Bn 

Mo = Mo = 2^1, Ml = m'i =0, S ^ Cn 



(5.50) 



^1 = ^0 = m'i = 0, S = B, 



MO = i'l , 

^0 = Mo = 2^1, Ml = m'i = 0, 



S = C, 



,"(5.51) 



(with A e P+). 

Let us now turn to the corresponding difference operators. Let Di^(r^s):- 



denote the operators Di,...,D„ (2.6) with parameters given by the table. 



We 



claim that the polynomials PA'.(ii,s), ^ £ 'Pji, are joint eigenfunctions of the operators 



,{R,S)- 



For R = (B) C'n, and for the polyno mials (5.5C), this is an imme- 



diate consequence of Theorem 3.1C 
By conjugating Di^(b„,s) with muj„ 
operator Di (2.10) with parameters 




For the polynomials (5.51) this is seen as follows. 
|5.47), one obtains (up to an additive constant) the 



= ^ Vb{eXj)Y\_Va{£Xj + Xk)Va{eXj 



- Xk 



l<j<n 
£ = ±1 



cosa{if3 + exj) -^pe 



^ - 1 



(5.52) 
{a = 1/2) 



Ecos a{ip + exj] 
cos aexj 

l<j<n 
£ = ±1 



Vb{exj 



) Y[ Va{exj + Xk)va{exj - Xk) ^ + const 



(calculate the residues to verify the second equality). Therefore, the polynomials ( p.5l[ ) 
are eigenfunctions of -Di,(_b„,s). One generalizes this to Dr,(B„,s)j r > 1, via similar 
reasoning as in the proof of Proposition 5.2: first one shows that D^ ^^^ 5) leaves invariant 



the space of anti-periodic polynomials mij^A (by calculating the residues); from the 
asymptotics for Im x 00 in th e positive Weyl chamber, it then follows that the operator 
is triangular. One uses Eq. (5.52), the monotony of the spectrum of D\.(b^ .gi (Le mma 5.1), 
and the commutativity of the operators to conclude that the polynomials (5.51) are joint 



eigenfunctions of -Di,(s„ 



, Dn,{Bn,S)- 



The value of the additive constant in the r.h.s. of (5.52) can be easily obtain ed b y 
spectru m of the operators on both sides of the equation (cf. eqs. (3.33)- 
1, and (p7|)): 

" (ch/3(M, + l/2)-ch/?pj) (5.53) 



const ■ 



l<3<n 



with pj = {n — j)a + ki/2 (S* = B„) or p. 
could generalize (5.52) to an expression tl 
Di, . . . , Dr with parameters (5.49). The p 
by co mparing the spectrum of the operato 
( |3.7l| ) in terms of £i,n(6'), . . .,Er,n{0) with 

Remarks i. The operator I?i,(_r,s) coincid 
Macdonald difference operator that is 
uji = ei of . For technical reasons, Macdi 
and the weight lattice of 5^ is scale d co i 
consequence that in comparing with |M2| 
factor 2 if S = Cn and = B(C)„. Specific 
and -Di,(fi,s) ~ -D2ei if 5 = C„ and R — B( 
ti. The operators I)i,{j5,5) , • • • , -C)„,(fl,s) 
difference operators that are simultaneous 

S = Cn and consider the operator 



...,e„ = ±l l<j<n 



with 



sin a{vi + 



sma{i'2 

This AAO coincides up to a multiplicative ■ 
is associated with the nth fundamental we 
One has (cf. Remark i) D'^ ^ ~ E^i^ if 

It follows from [M2] that D'^ ^-jj f^^-^ is di 

ever, that D'^ is not in P but its sqi 

eigenvalues and using Theorem 3.12). 

lii. From a group-theoretic perspective 
isomorphic to M['P]^ (the W-invariant pa 
r = Z"). If S = Bn, then one has Pgv = T 
P is a subgroup of Pgv with index two; so 
In the latter case, one can extend D to an a 
by replacing the generator Dn.(R,c„) by D[ 

5.4 The Root System D„ 

We conclude by briefly sketching the state 
should not have much difSculty to supply r 
subsection). Put 

MO = Ml = 

Then di,{z) — 1 and A(3;) (3.5C) reduces to 
the correspondence of parameters is via E( 
Macdonald's operator D^^^ associated with 
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For R = D„, the torus T r = M." / ) is the same as for R = _B„. The Weyl group, 
however, is smaUer: only an even number of sign flips of the variables Xj, j — 1, . . . ,n, is 
allowed. For k — 1, . . . ,n~2, the fundamental weights ujk of D„ are the same as in (5.42) , 
but ujn~i and are now given by the half-spin weights 



(£! + ■■■ + 



Cn-l 



e„)/2. 



(ei + ---+e„_i+e„)/2. 



(5.57) 



It is not hard to see that the cone of dominant weights V]^^ generated hy ut, k 
consists of the vectors 



= 1, 



(A + 5uj„), = (Ai + 5/2, A„_i + 5/2, e(A„ + 5/2)) 



with 



AgP+, 5 = 0,1, 
The Macdonald polynomials Px'^DnJ ^' G 



e = ±1. 



(5.58) 



(5.59) 



form an orthogonal basis of 



D„, Ao^dx). By combining the polynomials associated with (A + Slli„)+ and 
(A + S(jj„)- one obtains polynomials that are even in Xj, j = 1, . . . , n. These are related 



to Koornwinder's polynomials in the following way (cf. Eqs. (5.5C), (5.51)) 



P\+,D„ +P\_,D„ = Pa, 

P(X+ui„)+,D„ +P(A+^„)_,D„ = PA, 



MO = Ml = Mo = Ml = 0, 
Mo = Ml = Mo = m'i = 0, 
MO = Mo = m'i = 0, Ml = */3- 



if A„ = 0, 
if A„ > 0, 



(5.60) 

In the second and the third line of the above formula one obtains a sum of D„ polynomials 
rather than the polynomials themselves. Nevertheless, Eq. ( 5.6C ) determines the D„ poly- 
nomials uniquely. This is because flipping the sign of one of the a;j's in P(a+<5ij„)^,d„ results 
in P(A+(5t^„)- ,-Dn • Consequently, the coefBcients of the expansion of P(A+<5cj„)e,D„ in D„- 
typ e mo nomial symmetric functions are determined in terms of the coefficients occurring 

in (Eei). 



As regards the difference operators with parameters (5.56), the algebra D consists of 
commuting AAO's with the D„-type polynomials as joint eigenfunctions. One can extend 
D to an algebra that is isomorphic to 
Dn by 



E 



El- 



j=±l l<j<k<n 
= -1 



D' 



E 

ei,...,e„ = ±l l<:j<fc<n 



n 

:j<k< 

n 



jjvj'^'^" by replacing the generators Dn-i and 
va{e,x, + EfcXfc) e-«-i^i+-+-"^")/2, (5.61) 
vaie.x, + s,Xk) e-«^i'^i+-+^"^")/2. (5.52) 



These operators are proportional to Macd onald's operators Euj 
associated with the half-spin weights (p.57|). 



and , which are 



Appendix A: Cancellation of Poles 

In this appendix we prove two results, which were needed to dem onstrate that D,- maps 
into itself. It was claimed in the proof of Proposition |3.3| (Section p.2|) that the 



following expression: 



(-1)= n 



0gJlC...Cj^=j 
1<S<|./ 



(with A G P+ and J C {1, 



,n}, Jo = 0; 

and at xi — —Xj — 27, j — 2, . . . ,n (poles c 
to zeros in the denominators of the coeffici 
that we assume that the parameters 7, /i, 
are chosen in such a way that these poles 
the above regularity claims, thereby compli 
Before turning to the details, let us outlin 
of a sum of terms of the type 



where the index sets B^/ C J, s' = 1, . . . , s 



The terms (A.2) are associated with the sei 
C C C . . . C J 

(with the cell J fixed). Each term in (^ 
construct an involutive operation a (a^ =h 
a way that the terms associate d wit h a s( 
residue. Therefore, the poles in (A.l) cancf 



Lemma A.l (pole of type /) 

Let 7, /i, ^s, u's ( 5 = 0,1) and X2, ■ ■ ■ ,x„ bi 
poles. Then (A.l) is regular as a function 1 



Proof 

First note that the lemma is trivial if 1 ^ ^ 
xi. But if 1 £ J, then V} gives rise to a po 
Assume 1 G J and let Bs^ denote the bit 
define the following map a on the coUectioi 

l.A If \Bsi I > 1, then a maps (TO) to the 



0g Ji g J2 g ■■■ g Jsi-i 

l.B If \Bsi I = 1 and si > 1, then a maps I 
C 7, C /„ c . . . c 

2. If \Bsi \ = 1 and si = 1 (i.e. Ji = {1}) 
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Phrased in words: unless = {1}, the map a puUs the index 1 out of Bs^ and places 
it in a newly created block, which is sandwiched between Ba^ \ {1} and -Bsj+i (case l.A); 
when contains only the index 1, then a merges the b lock s Bsi = {1} and Bs^-i if 
si > 1 (case l.B) or, if si = 1, then it leaves the sequence (A. 4) unchange d (ca se 2.). 

Thus defined, a is indeed an involution on the collection of sequences (A.4): the cases 
l.A and l.B are inverse to each other (see Fig. 1. below). 



At xi = -7 the r.h.s. of ( |A.ll[ ) and ^AA 
xi = —7 cancel. 

If we are in situation 2., i.e. Bi — {1}, th 
in Vj is compensated by a zero in the differ 

[mx{x + 2761) - 



l.A 





jl,- ■ -jp, 1 









l.B 





jl, ■ --jp 


1 











This shows that the total residue at xi — 
the proof of the lemma. 



Fig. 1. A graphical representation of the map a. 

W e claim that in the first sit uatio n (i.e. l.A or l.B) the pole at xi — —7 in the term 
( |A.!^ ) (which is associated with (|A.4|) ) cancels against the pole in the term corresponding 
with the (T- image of the sequence ( |A.4| ) . To see this, we may assume th at we are in situation 
l.A. One obtains the term corresponding to the sequence (A. 5) from (A. 2) by making the 
substitutions 



mx{x + 2763 J 



s + 1, 



mx(x + 2'yeBi -27(5i,siei) 



(A.7) 
(A.8) 

(A.9) 
(A.IO) 



( denotes the Kronecker symbol). This substitution in (A.2) amounts to replacing the 
part 

{Vljvl^^^i}) [mx{x + 2^eB^)~mx{x)]^ (A.ll) 
^a{xj + xi) Va(xj + 3^1 + 27) [mx{x + 2jeBi) - mx{x)] 



jes.i\{i} 



by 



[mxix + 2jeBi - 275i,3iei) - mx{x)] = (A. 12) 

- Y\_ Va{xj + xi)va{xj ~ xi) [mx{x + 2^631 ~ 2'y5i^siei) ~ mx(x)]. 



Lemma A.2 (poles of type II) 

Let 7, ^, fis, ii'f, ( 5 = and X2, ■ ■ ■ ,x„ bi 
poles. Then (A.l) is regular as a function i 



Proof 

The proof is very similar to that of Lemma | 
if J does not contain the pair becau! 
xi = ~Xj — 27. 

Assume for the remaining part o f the p 
map on the collection of sequences {^-A)^. 



l.A If the pair {l,j} is contained in one 
\Bs^ I > 2, then aj maps sequence (A 



C 7, C . . . C 7 , C 7 



l.B If one of the bl ocks of sequence (A.4 
maps sequence (A.4) to 



g Jl g ■ ■ ■ C J, 

2. If _Bi = {l,j}, then Gj maps sequence 

3. If the pai r jl ,.?'} is not contained in 
sequence (A.4) to itself. 



.7SS=l\{l} 



It is clear that is an involution, the cases 
2. below). 
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l.A 









ji , ■ ■ ■ > 






l.B 


1,J 


-> 2. - 











Fig. 2. A graphical representation of the map aj. 

Consider situation 1., assuming case l.A. The a ppli cation of aj boils down to making 
the following substitutions in the associated term (A.2): 

(A.16) 

.} (A.17) 



s + 1 



■,J\Js 



These substitutions amount to the following change in the term (A.2): replace the part 

(A.19) 



by 



Va{Xk + Xl) Va{Xk + Xl + 2'y) Va{Xk + Xj) Va(Xk + Xj + 2'y) 
Y[ t^aixk + 3:^1) Va{Xk — Xl) Va{xk + X j) Va{Xk — Xj). 



(A.20) 



keBs ■ 



At Xl = —Xj — 27, (A.19) and (A.2C) differ o nly by sign. Consequently, the residues at 
Xl — —Xj — 27 of the corresponding terms in ( LV.in ad d up to zero. 

In situation 2. the pole in the coefficient of ( |A.2[ ), which is caused by Vg^, cancels 
against the zero in the difference of the monomial symmetric functions: 



[mx(x + 27e{i, 



J}) 



mx 



-27) 



= 0. 



In situation 3. the denominator of the coefficient of (A.2) has no zero at a;i 
so the term is regular at xi = —Xj — 27. 



(A.21) 

^j-27, 



We conclude from the above analysis that the total residue &t xi = —Xj — 2'y in the sum 



(AT) is zero, thus completing the proof of the lemma. 



subject to the condition 

1 < ji < j2 < • ■ • < is < n, 



1 < ii < 12 < 



< ir 



< n — r + 1. 



(B.8) 



We shall now show that these monomials cancel in pairs. 

Let a b e th e following operation defined on the above collection of pairs (B.7) with 
condition (B.^): 



A. If ii < ji or s = 0, then 

{(ji,j2, • • • ,js), (ii,i2, 

B. If ii > ji or s = r, then 

{{ji,h, ■ ■ - Js), {ii,i2, ■ 



,ir-s)} {(ii, ji 



,js), {i2, 



)}, 



,ir-s)} {(i2,...,js), 



(B.9) 



(B.IO) 



Roug hly speaking, a compares the first entries of the two elements constituting the pair 
(B.7) and moves the smallest of these two to the first entry of the other element. One 



easily verifies tha t: fi rst, a is well defined in the sense that the image of (B.7) is again 
a pair satisfying (B.8); second, a is an involution [a^ — id), the cases A. and B. being 



inverse to each other. 

For the associated monomial (B.6), acting with a amounts to an increase (case A.) or 
a decrease (case B.) of the number s by one, i.e. it fiips the sign of the corresponding 
monomial. Therefore, combining the term (B.6) associated with a pair (B.7) with the one 



associated with its image under a entails the vanishing of the sum (B.5), which completes 
the proof. 



Remark. If one replaces the upp er bound n — r+l of the second summation in (B.5) by n, 
then, for r = n, expression (IB. El ) also vanishes. (Indeed, the above proof again applies). In 
this case the vanishing of ( |B.5| ) amounts to a well-known relation bet ween the elementary 
symmetric functions and the complete symmetric functions (see e.g. [V[l|). 



Lemma B.2 The function 

Er ^71 (^1 5 ■ • ■ 5 in , Pr 7 ■ • ■ 5 Pn ) — 



0<3<r 



(-ir 




Pil ■ ■ -Prr-s 



(B.ll) 



1 < r < n. 



is the unique solution of the recursion relation 

Er^nitx, • • . ^ in, Pr ; 



Pn) = {tn — Pn)Er-l.n-l{tl, . . . , tn-\\Pr, ■ ■ ■ ,Pn) + 
Er,n-l{tl, . . . , t„-i;Pr, ■ ■ ■ ,Pn-l), 

l<r <n (B.12) 



with the convention 



Eo,n = 1, 



Er,„ = if n < r. 



(B.13) 



Proof 

It is clear that (B.li) with conditi on ( | B.1^ 
n). After splitting up t he sum in (B.ll) ir 
indeed solves Eq. (B.12): 



i. terms with js = n: 
li. terms with js < n and ir-s — n: 
Hi. terms with js < n and ir-s 



< n: 



Remark. In some cases Lemma B.2 can be 



For instance, one easily verifies with the aic 



Pr,Pr+l, 



then 



In particular, i5„,„(ti. 



{1 "} \jeJ 

\J\=r 

tn;Pn) = {tl - P, 



Appendix C: D{A^) = : 

In this appendix we present a result due t 
an AAO or PDO is zero on all sym met ric 
zero. This fact was used in Section B.5 to 



(Theorem 3.11), and again in Section \l.3\ 
BCn-type hypergeometric PDO's of Heckir 
Let 



,K1 ,K 



with K 



[Ki, . . . , Kn) m IK" or N'' m 
AAO's/PDO's of interest are of the form: 



D 



E 

l<r<J 

A' 



Vr-.UCR" ^ 

continuous on an open dense set U C M". 

Proposition C.l (D{A^') = ^ £» = 0) 
Let D be an AAO /PDO of the form (Iclj 



then 



Dmx = 0, 

Vl{x) = V2{X) ^ ■ ■ ■ - 
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Proof {\R3\) 



Introduce the following vector-valued functions: 

tx-.U^C^', {{f^a)mx){x),...,{f^iMymx){x)) , A G P+, (C.6) 



v.U 



^ (Vi{x),...,Vm{x)) 



(C.7) 



The fact that m\ is in the kernel of D translates itself geometrically in the orthogonality 
of tx and v: 

bmx = G^=^ {tx,v)=Q. (C.8) 
We will assume w ^ and derive a contradiction. Let A'^', . . . , A'*^' be vectors in ■ 



One has v _L t-^{i) , . . . ,t^{M) ■ Therefore, the vectors (x), . 



{x) must be linearly 



dependent for all a; € W for which v{x) 7^ 0. Since v{x) is continuous in x, there exists an 
open ball B <ZlA on which v{x) ^ 0. The fact that the vectors tx{s){x), 1 < s < M are 
real-analytic in x then entails 



det (Col[t;^(i) {x),..., t^iM) {x)]) = 0, 



yx£U. 



(C.9) 



We will now show that an appropriate choice of the vectors A'^', . . . , A'^^' contradicts the 
vanishing of the above determinant. 
Let A e P+ and y G M" be such that 



Ai > A2 > ■■ ■ > A„ > 0, 



yi > 2/2 > ■ ■ ■ > yn > 0. 



From the asymptotics (cf Eq. (B.37)) 

if^mx){x)U=^Ry^r^,xe^^'''^'', 



with 



one derives: 



lim e'-^'''^ tx{x)\x=iRy = (r^(i) 



R ^ 00 

(AAO) 
" (PDO) ' 



Pick the vector A (subject to condition (C.IO)) in such a way that 
■^kM.a / •^«(p),A' l<r<p<M. 



(C.IO) 

(C.ll) 
(C.12) 
(C.13) 

(C.14) 



That such a A exists follows in the AAO case from (C.12) and the fact that the vectors 



AM) 



are distinct; in the PDO case one can pick distinct prime numbers for the 



components of A. 

We use A = (Ai, . . . , A„) to form the vectors A'^', . . . , A'**^^ in the following way: 



A(^) = 



(s-l)(Ai,...,A„) 



((Ai 



,(A. 



On the one hand, Eqs. (CS) and (C13) imply 



■^k(i),a(i) 



(AAO) 
(PDO) 



^k(1),A(M) 



1< s < M. 



^k(m),A(1) '^k(a^),A(m) 



(C.15) 



(C.16) 



On the other hand, for the above choice ol 
determinant: 

Therefore, 

- n < 

l<r<p<A/ 

Because A is chosen such that r,.(,.i n 7^ t.,(„ 
determinant r 0, contradicting ( C.16 ). 
Hence, v (C.7) must be zero. 
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